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8 ■ Abstract 

Dirac-Born-Infeld type effective actions reproduce many aspects of string theory classical tachyon 

C$ ; 

dynamics of unstable Dp-branes. The inhomogeneous tachyon field rolling from the top of its 
potential forms topological defects of lower codimensions. In between them, as we show, the tachyon 
energy density fragments into a p-dimensional web-like high density network evolving with time. 
We present an analytic asymptotic series solution of the non-linear equations for the inhomogeneous 

O ■ 

■ tachyon and its stress energy. The generic solution for a tachyon field with a runaway potential in 

o 

"sj- I arbitrary dimensions is described by the free streaming of noninteracting massive particles whose 



initial velocities are defined by the gradients of the initial tachyon profile. Thus, relativistic particle 
mechanics is a dual picture of the tachyon field effective action. Implications of this picture for 



inflationary models with a decaying tachyon field are discussed. 
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Introduction. In this paper we investigate generic inhomogeneous solutions of Dirac-Born- 
Infeld type theories 

S = - J d p+1 x V{T)\Jl + a'd^Td^T + O (d^T) , (1) 

where T(x^) is a (dimensionless) scalar field and V(T) is its runaway potential (no minima). 
In string theory a' is a square of the fundamental length scale; we put a' — 1. The action (JT|) 
was proposed in |l| as an effective field theory description of the open string theory tachyon 
which describes unstable non-BPS D-branes. In application to the string theory tachyon (JT|) 
should be understood in the truncated approximation, i.e. valid only in the regime where 
higher derivatives are not large. The potential is often chosen to be V(T) = r p /coshT 
for the bosonic case which we consider. At large T the potential has a runaway character 
V(T) ~ e~ T with the ground state at infinity. 

There are several motivations for studying the properties of the effective action It 
is difficult to find the open string tachyon dynamics for generic tachyon inhomogeneities. 
The action meanwhile, permits us to study complicated tachyon dynamics in terms of 
classical field theory. The relatively simple formulation of tachyon dynamics in terms of the 
effective action (P) has therefore triggered significant interest in the investigation of the field 
theory of the tachyon and the possible role of tachyons in cosmology. Indeed, the end point 
of string theory brane inflation is annihilation of D — D branes, which leads to the formation 
and subsequent fragmentation of a tachyon condensate Q]. Thus the potential role of the 
tachyon in cosmology cannot be understood without first understanding its fragmentation. 

Apart from its application to the string theory tachyon, the search for the structure 
of general solutions of the theory is an interesting mathematical problem in and of 
itself. The equation of motion arising from the action (0) is an example of a complicated, 
non-linear, partial differential equation which, as we will show, admits a relatively simple, 
general, inhomogeneous solution. The evolution of the tachyon field T(t, x) can be viewed 
as a mapping T(t ,x ) — > T(t,x) that becomes multivalued and generates singularities at 
caustics |3|. Besides the DBI type theories (JTJ), there are other cosmologically motivated 
phenomenological models of fields with high derivatives which share the problematics of (0). 

If the potential V(T) is symmetric around T = and the inhomogeneous tachyon field 
begins rolling from the top of its effective potential, then topological defects (kinks) can 
form due to symmetry breaking 4] . In this paper we consider what happens to the tachyon 
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FIG. 1: The focusing of energy density into a web-like structure due to caustic formation. Whiter 
regions correspond to higher densities. The left panel shows a nearly homogeneous initial Gaussian 
random field profile and the right panel shows the same field a short time later. 



field in the region where it rolls down one side of the potential. We will see the formation 
of sharp features in the tachyon energy density due to its fragmentation. These features, 
which are related to the convergence of characteristics of the field T, have to be distinguished 
from topological defects. The full picture must incorporate both effects, formation of kinks 
and tachyon fragmentation in the space between them. Because it is hard to study with 
CFT tachyon dynamics with a generic, spatially varying profile |5|, previous calculations 
dealt with a plane wave tachyon profile [6] . In this case the tachyon decays into equidistant 
plane-parallel singular hypersurfaces of co-dimension one, which were interpreted as kinks. 
The effective action for a plane wave tachyon predicts similar result, as we will see later. 
However, this inhomogeneous profile is atypical in the sense that fragmentation between 
kinks does not occur. In the general case we expect both structures, weblike fragmentation 
and topological defects. 

In this paper we concentrate on tachyon fragmentation between kinks. We begin by 
showing an image that illustrates the fragmentation of the tachyon field as it rolls down one 
side of its potential. Figure^shows the result of a numerical lattice simulation of the energy 
density of a two dimensional tachyon field rolling down one side of its potential, as described 
by the equation of motion (J2J) below. We used the LATTICEEASY code pfl adapted for eq. 
(J2J). Starting from an initial random Gaussian field T the energy rapidly became fragmented 
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into an anisotropic structure of clumps joined by filaments into a web-like network. 

The tachyon energy density pattern in Fig. ^ is reminiscent of the illumination pattern 
at the bottom of the swimming pool or the web-like large scale structure of the universe. 
The similarity is not coincidental: the underlying mathematics has common features in all 
three cases. In what follows we present an analytic solution to the tachyon equation of 
motion that describes in detail the formation of this structure. We begin by describing 
a good approximation to the dynamics of eq. (J2J) and go on to show how to extend this 
approximation into an asymptotic series for the field T. 

The Free Streaming Approximation. The equation of motion for the tachyon field follows 
from the action (|T|) 

d u d»T d „ TgrT _ Kt = Q , 2 s 

" l + d a Td a T V v ; 

For simplicity we will confine ourselves to a pure exponential potential V(T) = e~ T , however 
our results are qualitatively valid for any runaway potentials. The energy density of the 
tachyon field p = T 00 is 

-T 

e 



P = i T + e- 1 J I + d,Td»T , (3) 

j l + dpTdxT V 



We have observed solving (j2J) numerically jsj that if we define an operator P(T) = 
1 + dpTd^T the field T rapidly approaches a regime in which P(T) P3 0. We write this by 
saying 

T(x^) ss S(x») , (4) 

where S satisfies the equation 

s 2 -(y x sf = i. (5) 

The dot represents a time derivative and the spatial derivatives are with respect to the p 
spatial coordinates x on the brane. This equation is the Hamilton- Jacobi equation for the 
evolution of the wave front function of free steaming massive relativistic particles. 

Let us consider this particle description. At some initial time to we can label the position 
of each particle with a vector q. Equivalently we can say that q parametrizes the different 
particles. The initial four-velocity of the particle is given by O^Sq. If we further define 
the proper time r along each particle's trajectory, we can switch from coordinates (t, x) to 
(r(i, x), q(t, x)) and obtain an exact parametric solution to (jSJ 

x = q - VjSqt (6) 



t = yfl + \VjS \ 2 T (7) 

S = S + r. (8) 

The interpretation of the solution © is very simple and intuitive. It tells us that the field S 
propagates along the treajectories of the massive relativistic particles, growing linearly with 
proper time. The slope of each characteristic depends only on the initial gradients of So on 
that characteristic. 

In geometrical optics photons are massless, but the qualitative picture of their wave front 
propagation is similar. This expains the similarity between the two dimensional web-like 
pattern of Fig. ^ and the illumination pattern at the bottom of a swimming pool. The 
focusing of particle trajectories corresponds to higher density concentrations and further, to 
the formation of caustics at the loci where trajectories cross. 

The Full Solution. Let us now consider the energy density. Looking at equation (j3J) we 
see that the exponential pieces are growing exponentially small, as are the arguments of 
the square roots. The second term in Q will thus rapidly become irrelevant and we need 
consider only how the exponentially small numerator and denominator of the first term will 
be related. (The leading term in T 2 is one). 

To calculate p we will need to go beyond the free streaming approximation (j3J. In view 
of the exponential in the numerator of p we conjecture that T can be expanded as 

T^^S + MSJe" 25 , (9) 

where fi is a sub-exponential functional of S. (We could include a lower order term foe~ s . 
As we explain below, we can solve exactly for f in that case and we find that its effects can 
be absorbed into S and j\.) Now we are going to check the validity of this expansion. 

Plugging the expansion (J2j) into equation (j2j) and keeping only terms proportional to e~ 2S 
gives the following equation for the sub-exponential function f\ 

(d»Sd u S) d^h + 2 (1 - US) & Sd„h ~ 4nS A = , (10) 

where US = d y d u S = -S + V/S. 

This equation can be dramatically simplified by changing from (t, x) coordinates to (r, 
q) coordinates, as defined by the characteristics of S in (j^HZj)- In these coordinates, fi has 
no spatial derivatives and equation (|10|) reduces to 

/ 1)TT - 2 (1 - US) f 1)T - 4DS h = 0, (11) 
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Note that OS can be calculated either with respect to (t, x) or (r, q) coordinates. We can 
further simplify this equation by introducing a new variable y 



y = 2f 1 -f ljT . (12) 

Equation (jllj) can be rewritten in terms of y 

y, T + 2OSy = 0, (13) 
which can be immediately solved to give y(r, q) = yo{q) exp (—2 J T dr' OS). From this and 

fi(r,q) = fu(q)e^ f dr' e~ 2 ^ 2 f ' ** ' ns . (14) 

To proceed further we need to calculate OS. In principle it can be done from the solution 
El) in parametric form by inverting the (t, x) — > (t, q) coordinate transformation matrix 
12j. Instead we will use the following trick. Plugging the expansion Q into the energy 



density Q we find p ~ 1/y 2 (2/i + d^Sd^fi). Observe that the denominator here is exactly 
identical to y/2y, thus 

P ^ ^ = Po(q) ef . (15) 



This is precisely the expression for the energy density of free streaming relativistic particles 
which obey the relativistic continuity equation d^S (pd^S) = 0. In fact, ([To]) is the solution 
of this continuity equation in the coordinates {r,q). However, there is another form of the 
energy density, which is equivalent to (fT3|) 

\dg\ 

where the denominator is the Jacobian of the x — > q transformation. Indeed, from conser- 
vation of the energy density in a differential volume we have d p qpo(q) = d p xp, which leads 
to the formula (|16|). It is now straightforward to calculate the Jacobian from the formulas 
(JBHH). In p dimensions it is a polynomial in r of order p. For example, in the one dimensional 
case |||| = (t c — r), where r c (q) is a function of the gradients of So. In the two dimensional 
case |||| = (t c1 — t)(t c2 — r), where r cn are function of Vg-So- 

Comparing expressions ([To]) and (|TB|) . we can find OS and calculate f\ using (JT3J). We find 
that a constant originating from the integration in |T4|) can be absorbed in S while the rest 
of fi will be a 2p order polynomial in r. This demonstrates the validity of the approximation 



(JHJ) by showing that the f\ term provides only exponentially suppressed corrections to the 
leading term. 

We could go further and include other powers of e~ s in our expansion (jUJ). 

oo 

T(xn^S + Y,fn[S]e-^ s . (17) 

n=0 

We have explicitly checked all such terms up through e -45 , including a possible term propor- 
tional to e~ s , and found that they simply provide corrections to the integration constants 
of S and f\ plus terms that are exponentially suppressed relative to the ones we have dis- 
cussed. We can further show that all such terms f n have the same characteristics as S, 
and we therefore conclude that the general inhomogeneous solution T propagates along the 
characteristics (07J). Up to exponentially small corrections (which could in principle be cal- 
culated order by order), the complete solution for T can be represented by the two functions 
S (q) and f u {q). 

With these results, we can evaluate the energy density ([Tfi|) to leading order using only 
/i. For an arbitrary brane dimension p we have 

P^PO Ua-Xnty 1 , (18) 
n=l 

where X n (q) are the eigenvalues of <9 M (S^d^So^. From here we see that the energy density 
first reaches large values in regions where X n (q) is maximal. For some critical trajectories 
q c at a critical time t c the energy density becomes singular, which corresponds to caustic 
formation. This is exactly what one would expect in the picture of free streaming, massive, 
relativistic particles, where the energy density blows up at the orbit crossings. 

Conclusions. Our most important conclusion is that the general inhomogeneous solution 
of the field theory (JIJ) very rapidly approaches the asymptotic form (|17|h which is equivalent 
to the relativistic mechanics of freely propagating massive particles with velocities v a = 
— V fSo(q a ) ■ In other words, there is a duality between the two Lagrangians 

V(T)y/l + d»Td»T J2 V 7 * 3 ^ " ( 19 ) 

a 

The whole process of unstable Dp brane decay in the the dual picture is described as "crum- 
bling to dust" of massive particles. A complete discussion of the interpretation of massive 
particles and anisotropic high density structures (clumps, filaments, sheets) which they form 
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would be beyond the scope of this paper. In the 1 + 1 case the high density regions of the 
orbit crossing were conjectured to be DO branes 

From the velocities of the characteristics (JBJ) we can see that around maxima of the 
initial field profile the characteristics will tend to diverge and the profile will flatten. In re- 
gions around minima, however, characteristics will tend to converge, the profile will become 
sharper, and after some critical time t c the field solution will become multivalued. In the 
dual picture of relativistic particles this corresponds to caustic formation. At the caustics 
the energy density blows up. Caustic formation also entails divergences in the second deriva- 
tives of T, which signal the breakdown of the truncated approximation (0). In short the 
Lagrangians (|19|) are unable to describe the field T when its solution becomes multi-valued. 
In the picture of freely moving massive particles we can inlcude interactions to cause them 
to stick together as their trajectories intersect with impact parameter ~ \fa! '. 

Let us make a remark about a one-dimensional plane wave tachyon profile T = cos(a;) 
rolling from the top of its (symmetric) effective potential. Since the parts of the field that 
roll to the right have no minima they do not form caustics, and by symmetry the parts 
rolling to the left do not either. In this particular case the tachyon fragments into kinks 
only. 

We can also consider a more general profile, however. Tachyonic instability occurs for all 
inhomogeneous modes k for which the effective mass m 2 = k 2 — 1/a' 2 is negative. Therefore 
the generic tachyon initial profile is a superposition of a number of modes, which produces 
a random Gaussian field T (q). These initial conditions typically arise from quantum fluc- 
tuations during symmetry breaking (see e.g. j^]). In this case we once again expect the 
formation of topological defects due to symmetry breaking. Outside these defects, however, 
the tachyon field will fragment into a web like structure as shown in Fig. ^ If the defects 
are walls these webs will form within each domain; if they are strings the web of caustics 
will be mixed in with the strings. 

If the spectrum of T (q) inhomogeneities has scaling properties (as quantum fluctuations 
do), then the web-like network will evolve in a scaling manner. The smallest scale of the 
web is defined by the largest tachyonic mode k. The dual picture of freely moving massive 
particles which stick together as their orbits intersect gives a simple explanation of such 
fragmentation. 

Finally, we note the relevance of our result for cosmological applications of the tachyon. 
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In the context of brane inflation, which ends with a pair of D3 — D3 branes annihilating, 
the tachyon is a complex field and strings will be created and the rst of energy is transform 
into radiation Q|. If the result of real tachyon field fragmentation which we derive above 
is extended to the complex tachyon, then annihilation of D3 — D3 branes results in the net 
of strings plus massive particles with the matter dominant equation of state. The absence 
of radiation domination after brane inflation may pose a problem for the model 
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